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Abstract 

We propose a very simple but general method to construct solutions of Maxwell 
equations propagating with a group velocity v gr / c. Applications to wave guides 
and a possible description of the known experimental evidences on photonic tunnel- 
ing are discussed. 



1 Introduction 



It is a well known tool that solutions of Maxwell equations in vacuum 

V • E = VxE = 



V ■ B = V x B 

are electromagnetic waves: 



1 dB 

~c~dt 

(1) 

1 <9E 

c ~dt 



4> being any component of E,B. From eq. (Q) follows that the velocity of advancement 
of the related wave surface (front velocity) is c [[TJ], the speed of light in vacuum. For a 
wave packet, it is of particular relevance also the concept of "group velocity", being the 
velocity v gr with which the maximum of the wave packet propagates. To this regard, it is 
commonly believed that the group velocity of an electromagnetic wave in vacuum takes 
always the value c. However, there is no proof @ that electromagnetic waves with v gr = c 
are the only solutions of eq. fll]) or (|2]) and, in particular, there are many works || |]] 
dealing with both v gr < c and v gr > c solutions. 

On the experimental side, many recent experiments on photonic tunneling conducted 
with different techniques and in different ranges of frequency H || 0, f| have shown that, 
in peculiar conditions, (evanescent) electromagnetic waves travel a barrier with (group) 
velocity v gr > c. In these experiments, as remembered in |10[ , only the group velocity can 
be determined, so that these results obviously do not violate Einstein causality because, 



according to Sommerfeld and Brillouin [11], it is the front velocity (not group velocity) to 
be relevant for this and, as stated above, the Maxwell theory predicts that electromagnetic 
waves in vacuum have always a constant front velocity equal to c. 

From the theoretical point of view, the difficulty in the interpretation of the experi- 
mental results lies mainly in the fact that in the barrier traversal no group velocity can be 
defined, the wave number being imaginary (evanescent waves), so that the time required 
for the traversal (directly measured) is not univocally defineable pf|T2|. 



In this paper we want to give a simple but general method for building up solutions 
of Maxwell equations propagating with group velocity v gr ^ c that generalizes previous 
calculations |4[] and to show how the theory could describe, almost qualitatively, the ex- 
perimental evidence for v gr > c electromagnetic propagation. To this purpose, in the 
following section the general formalism is outlined while in section 3 this formalism is 
applied to electromagnetic signals in a wave guide. In section 4 the experimental evi- 
dence on photonic tunneling is reviewed and discussed, and its possible description in our 
approach is shown together with future possible tests. 



1 Incidentally, we quote the introduction of the concept of "phase time" which is a generalization of 
that of group velocity, and it is applicable in this case. 
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However, it is worth noting that the theory developed in section 2 can be applied as 
well to propagation of usual electromagnetic waves in a medium (making the substitution 
c — > c/n, n being the refractive index); in this case our approach can be an alternative 
one to well known phenomena such as, for example, plasma oscillations, Cherenkov effect 
and so on. 



2 General formalism 

In the experiments performed on superluminal barrier traversal, the wavelength A of the 
incident electromagnetic signal is always greater than or, at least, of the some order of 
a typical length of the experimental apparatus ||, [7], |[ ||. So, we have to use an 
approximation just opposite to the eikonal approximation, and we search for solution of 
eq.(||) in the form 

0(x,t) = M^t;k,oo)e^-^ (3) 

with an amplitude O strongly dependent on 4-position x, t. Substituting in eq.(^j) we find 
the equation satisfied by 0o : 

0Vo + 2i(^ + k.V0„) + (j!-^„ = O (4) 

For illustrative purpose we consider one-dimensional problems; taking O dependent only 
on £ = z — v t (and oo,k), where v is a given parameter (as well as u,k) that we will 
show to be identified with the group velocity, equation (^) now assumes the form 

/ v 2 \ d 2 (j) / vuj\d(f) (u 2 2 \ 

Here we observe that, if u = ck (and then v — c), eq.(^) is automatically satisfied, 
whatever being the dependence of 0o on 4-position; however, for O C-d e P en dent, these 
solutions are not the only ones. As we will show now, there are also v ^ c solutions; 
these can be obtained imposing a "quantification rule" (in the sense that it quantifies the 
effective dispersion relation) specifying how strong the eikonal approximation is violated. 
Let us impose that 0o in addition to eq. (|5|) satisfies also 

with a given constant Q 2 > 0. In other words, we costrain O to satisfy a harmonic motion 
equation whose frequency 7^7 (7 = 1/ a/1 — v 2 /c 2 ) quantifies the effective dispersion 
relation. In fact, substituting eq. (H) in @ we find that 

0o = 0' o e ±7fX (7) 

is solution of eq. @ with 

u 2 = c 2 (k 2 + n 2 ) (8) 
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and 

c 2 k I c 2 n 2 

< c (9) 



to V oj 2 



For real v we see that, for fixed lu, is a product of waves; in this sense, provided 
jQv <C oj and noted that 

(10) 



\U3 1 



dk 

v may be interpreted as a group velocity. 

Notice that for very small u and finite Q both k and v are imaginary, so that in this 
case, in a common language, we are describing evanescent waves. 

Now, let us consider a wave packet with a given spectrum in u 



[z,t) = I du<f/ (u)e** v)niz ^ iu)t) e i M u) '-^ ) = 

dw^(w)e i(4 ' z - w ' !) (11) 



with k'{uj) = k(u) ±7^2, uj'{u) = u ±jQv(u). Suppose that in the spectrum <p' (uj) of the 
wave packet one has a pronounced maximum ujq\ expanding the integrand in u around 
ujq we get 



■oz-uot) A^^^ e i{[± 7 n+|^| wo (a)-£Jo)]«-[±7n||| w + (w-a;o)]t} /j_ 2 ) 



i(feo 



(fc(a>o) = fco). From this, it is easy to obtain the group velocity of the wave packet: 



Vgr = Ihn dk |fcr r ~T = ^rUo = w (^o) (13) 



Let us stress that for fl ^ we have f (c^o) < c (see eq. @); we then obtain solution of 
Maxwell equations propagating with a group velocity lower than c depending on "con- 
struction" boundary condition (see below) through Q 2 . In fact, the space-time evolution 
of wave amplitude, and then Q 2 in eq. (|6|), can be "constructed" experimentally and the 
dispersion relation (§) or the group velocity in eq. @ can be further measured. 
Now, instead of eq. (|]), let us impose the following equation 

again with Q 2 > 0. Substituting in eq. @ we now easily find that 0o is solution of eq. 
([U]) with a dispersion relation 
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c 2 (A; 2 - fi 2 ) (15) 



propagating with a superluminal group velocity 



v = 



c 



CO 



■ 2 k 




■m 2 



> c 



(16) 



Note that even in this case <po satisfies a harmonic motion equation fll~4]) ; now, however, 
both k and v are always real, so that we have no evanescent waves. 

Incidentally, let us observe || that our subluminal solutions can be equivalently con- 
structed by requiring that <f) Q satisfies an Helmoltz wave equation on which solutions a 
Lorentz boost, for example in the z direction, is applied. However, we again point out the 
fact that a quantification rule can be realized only by construction of a given experiment. 

Let us now discuss this last point and, in particular, the introduction of the param- 
eter Q in the effective dispersion relations. For simplicity, we confine ourselves only to 
subluminal group velocities, but the same will remain valid also for the superluminal case. 

Our method is based mainly on eq. (||) which admits, for given initial conditions, an 
univocal solution (f> in (|7|) and then in This wave, product of waves, is what one can 
measure and is univocally determined by the experimental setup employed. In this sense, 
the parameter Q is given "by construction" once the experiment is given. For example, 
in experiments operating in the eikonal approximation regime (this is not the case of 
photonic tunneling experiments), namely when the wavelength A of the electromagnetic 
signal is much lower than a typical length of the experimental apparatus, the parameter 
Q is zero. On the contrary, in experiments with wave guides one can experimentally 
determine Q by measuring the (effective) "cutoff" frequency (as we will show in the next 
section, the introduction of Q leads to a modification of the cutoff frequency of the wave 
guide) . 

Finally, we stress that the actual theory can be consistently tested only in problems 
in more than one dimension. In fact, let us consider, for example, an incident signal 
on a wave guide: it is well known that the signal effectively propagates (with a finite 
group velocity) only in the axial direction, while across the guide a non-physical signal 
propagates with an arbitrarily high group velocity. In the following section we will show 
how our approach can be applied to realistic problems by giving a specific example. 

3 An application: wave packets in a wave guide 

As an application let us consider the propagation of an electromagnetic wave packet in 
a hollow wave guide, placed along the z axis, of arbitrary (but constant) cross-sectional 
shape with boundary surfaces being perfect conductors (the development of this paragraph 
is a generalization of chapter 8 of [[H|]). If uj is the frequency of the incident signal, let us 
write the spatial and temporal dependence of the electric and magnetic field inside the 
guide as 




(17) 
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where k is an as yet unknown wave number and ( is given in the previous section. The 
fields E ,B satisfy the wave equation 



(vi + (£-*.) + »(*-5)» + (i-*)*)(*)-o m 

where Vi = £ + gr- 

Now we impose the ^-dependence of the fields to be given by one of the quantification 



rules eq. (0) or (0); then the x, y dependence is determined by the equation 



+ ^ - * 2 t n 2 ) ( ! = o (20) 



c 2 



Bo 



where the upper sign refers to propagation with v < c and the other one to v > c. 
From this follows that, considering for example subluminal group velocities, in terms of 
the transverse and parallel (respect to the z guide axis) components, Maxwell equations 
become Q 

(k + jQ) E ± + - (uj + -fQv) e 3 x B ± = -iV ± E z (21) 

(k + 7O) B ± (u + jttv) e 3 x E± = -iV±B z (22) 

c 

e 3 ■ (V± xEi) = - {uj + -fttv) B z (23) 
c 

e 3 ■ (V ± xBj = - - (u + jQv) E z (24) 

c 

V ± • E x = -i(k + 7O) E z (25) 

V±-B ± = -i(k + 7fi) B e (26) 

e 3 being a unit vector in the z direction (for v > cit suffices to replace 7 = (1 — f 2 /c 2 ) -1 / 2 
with 7 = [v 2 / c 2 — l) -1 / 2 ). We deduce that if E Z ,B Z are known, then from the first 
two equations (pTD,(p2|) the transverse components of E,B are determined, once the in- 
dependence is known from the quantification rules; this leads, with respect to the normal 
case in which v — c, to an "effective" propagation frequency (and wave number) inside 
the guide given by uj' = uj + 7^t> (or oj' = uj + ^VLv for v > c), quantified just by the 
quantification rules. 

In fact, for TM waves, specified by the conditions 

B z = everywhere 

(27) 

E z = on the guide surface 



2 In general what follows remains valid if the substitution — >■ — is performed. For clearness of 
notation we restrict only to one sign. 
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we have 

E ± = ^V x i? 2 (2f 
while for TE waves, specified by the conditions 

i£g = everywhere 



(29) 



= on the guide surface 



dn 

the transverse components of the magnetic field is 



Bi = • V ± B z (30) 



In every case 
with 



B ± = Ye 3 x E ± (31) 



Y = — — TM waves 

c(k + 7^) 



c(fc + 7g 
r = — 1 E waves 



(32) 



In eqs. 
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r = — — A; =f (33) 



c- 5 

are the eigenvalues of the wave equation ( p0|) with the boundary condition given respec- 
tively by eq. (p7|) , (|29|) . For example, for a rectangular wave guide with cross-section a ■ b, 
the eigenvalues for TE waves are |13| : 

TL = rr 2 + (34) 



a 2 b 2 r 

with m, n positive integers. 

The equation ( p3|) is the dispersion relation for k: for an incident wave with a given 
u, the propagating modes are those with the wave number 

c 2 k 2 = uj 2 - c 2 (r 2 ± n 2 ) (35) 

First, let us consider the case in which v < c; the propagation inside the guide effectively 
happens only for real values of k, i.e. for frequencies uj greater than the "cutoff" frequency 



UJ C . 



cVr 2 + tt 2 (36) 



6 



while for u < lu c we are in the evanescent waves regime. Note that the effective cutoff 
frequency uj c is greater than the usual value u 0c = cT (fl = 0) and furthermore, for 
evanescent modes, the z-dependence does not show a pure exponential decay but is of the 
type e % ^ z e~\ k \ z (obviously we take Q always real). 

Now, let us consider the most interesting case for which v > c; we can easily prove 
that no cutoff arise. In fact, calculating the group velocity from eq. (|35|) 



u; 2 c c 2 Q 2 



c\/l--f + — (37) 
in 1 



we find that it is effectively v > c only if the condition 

c 2 tt 2 > cu 2 c (38) 

is fulfilled. This condition makes that in eq. (|35|) the wave number is always real. Note 
that if c 2 Q 2 < uj 2 c we again have propagation with v < c, with the same features as 
discussed above, but now the effective cutoff frequency u c = c y/T 2 — Q 2 is lower than 
the usual value u> 0c . 

The obtained results bring to a very important consequence: a group velocity (and 
then a traversal time) can always be coherently defined for v > c propagation. Notice 
that the fact that Q is a property of the experimental setup (see at the end of the previous 
section) guarantees the univocity of the identification v gr = v(u ) (v depending on Q) and 
so ambiguities in the definition of the traversal time are really not present. 



4 Discussion on photonic tunneling 



Recently Martin and Landauer |L4] have shown that the propagation of electromagnetic 
evanescent waves in a wave guide can be viewed as a photonic tunneling process through a 
barrier; this facilitates the experimental study of tunneling phenomena, due to the charge 
neutrality of photons with respect to other particles such as electrons. From the theoretical 
point of view, an open question is that of the barrier traversal time, because in the barrier 
region the momentum of the tunneling particle is imaginary so that no velocity can be 
defined. There are several approaches to the problem [|12j leading to different definitions 
of the traversal time. Instead, experimentally this time can be univocally measured, for 
example, from the coincidence of two photons, one travelling through the barrier and the 
other travelling in vacuum. 

Enders and Nimtz || ^] have studied photonic tunneling by means of microwave trans- 
mission through undersized wave guides operating below their cutoff (of the order of 6 4- 9 
GHz) and have obtained traversal times (for opaque barriers) from pulsed measurements 
in the time domain or, indirectly, in the frequency domain. Alternatively, Steinberg, 
Kwiat and Chiao |7j and later Spielmann, Szipocs and Krausz employed some ID pho- 
tonic band-gap material as barrier for measuring tunneling times, in the UV and optical 
region respectively, with the aid of a two-photon interferometer. 
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The experimental evidence can be summarized as follows. Both in microwave and 
in photonic band-gap experiments photonic tunneling is observed; if the barrier medium 
is non-dissipative, the traversal time in this opaque region is nearly independent of the 
barrier thickness (Hartman effect fllSl), so that for particular values of this length super- 
luminal group velocities have been inferred. Furthermore, two strange properties have 
been detected. First, in microwave experiments, dissipative tunneling studies have shown 
that the Hartman effect disappears with increasing dissipation ||. Second, the measure- 
ments of the tunneling of optical pulses through photonic band-gaps reveal that the pulses 
transmitted through a particular sample are significatively shorter than the incident ones 
||; this effect disappears for increasing transmission coefficients. The latter two effects 
seem to indicate, in our opinion, that real propagation in the opaque region happens, and 
this stimulates to apply our theory to the present case. For example, we stress that a 
usual mass term (for v < c) or a tachionic mass [16] for the photon (for v > c) cannot 
take into account the observed superluminal tunneling, because of the dependence on the 
barrier thickness of the measured group velocity (i.e. the dependence on a characteristic 
parameter and not on an intrinsic one such as usual or tachionic mass). Instead, our quan- 
tification frequency Q is not an intrinsic property, but would just depend on the employed 
experimental setup, so that the dependence on it of the group velocity (see for example 
eq. ( |37D ) seems to go in the right direction for taking into account the superluminal 
tunneling. Moreover, the simple fact that in our approach a group velocity can always be 
defined for superluminal propagation eliminates the ambiguities in the definition of the 
traversal time. 

On the other hand, our theory can be tested independently from photonic tunneling, 
for example constructing "ad hoc" an electromagnetic apparatus which realizes eq. @ or 
(|TJj) in some regions and then measuring the dispersion relation between u and k or the 
group velocity of the propagating waves. 



5 Conclusions 

We have studied the propagation with group velocity v gr ^ c of the solutions of Maxwell 
equations (in vacuum) and have shown that it is possible (and not violating Einstein 
causality) provided a peculiar space-temporal dependence of the wave amplitude is given 
through a definite quantification rule. 

As an application of the presented formalism we have considered the propagation in a 
wave guide and obtained, for v gr < c, a different effective cutoff frequency respect to the 
normal case, while for v gr > c no effective cutoff arises so that the wave number is always 
real and a group velocity can always be defined. In the evanescent regime for the v gr < c 
case, moreover, the real propagating waves show not a pure exponential decay but are 
only damped waves. 

Discussing photonic tunneling, we have pointed out how the presented approach can 
qualitatively describe the experimental evidences on this effect, even if definitive con- 
clusions are not yet reached and further experimental and theoretical investigations are 
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needed. 
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